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I. INTRODUCTION 

For quantum computers to become a reality we need 
to find or build physical systems that faithfully imple- 
ment the quantum gates used in the algorithms of quan- 
tum computation. The basic requirement is that the ex- 
perimenter has access to two states of a quantum sys- 
tem that can be effectively decoupled from environmen- 
tal noise for a sufficiently long time, and that transitions 
between these two states can be controlled to simulate 
a number of elementary quantum gates (unitary trans- 
formations). Systems that have been investigated inten- 
sively are atomic levels in ion traps [ll, H, superconduct- 
ing device physics using Joscphson rings [3| , nuclear spins 
4 1 (using NMR in suitable molecules) and quantum dots 
5|. In this paper we demonstrate the implementation 
of quantum gates using one-dimensional spin-J systems. 
The results are obtained using a computer simulation of 
these systems. 

The Hamiltonian of the XXZ model with kink bound- 
ary conditions is given by 



L-l 



H 



a=-L+l 



(1) 



where and are the spin-J matrices acting on 
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the site a. Apart from the magnitude of the spins, J, the 
main parameter of the model is the anisotropy A > 1 
and the limit A ^ cx) is referred to as the Ising limit. 
In the case of J = 1/2 kink boundary conditions were 
first introduced in Q. They lead to ground states with 
a domain wall between down spins on the left portion of 
the chain and up spins on the right. The third compo- 
nent of the magnetization, M, is conserved, and there 
is exactly one ground state for each value of M. Dif- 
ferent values of M correspond to different positions of 
the domain walls, which in one dimension are sometimes 
referred to as kinks. In Q, Koma, Nachtergaele, and 
Starr showed that there is a spectral gap above each of 
the ground states in this model for all values of J. Re- 
cently Q it was shown that for spin values J > | and 
for sufficiently large value of the anisotropy A the low 
lying spectrum of ([T]) for each value of M has isolated 
eigenvalues that persist in the thermodynamic limit. 

The presence of isolated eigenvalues is ideal from the 
point of view of quantum computation. The idea is to 
use the subspace, denoted by T>, of the ground state and 
the first excited state of the Hamiltonian to encode a 
qubit. We let the system evolve under its own unitary 
time evolution generated by the Hamiltonian ([1]) with 
the addition of a few local control fields. We have two 
requirements to fulfill: the time evolution should leave 
the qubit space P approximately invariant, and the (ap- 
proximately) unitary matrix describing the dynamics re- 
stricted to T> and stopped at a suitable time should co- 
incide with the desired quantum gate. 

The control inputs needed to drive the system such 
that high fidelity gates are obtained are determined using 
techniques from optimal control theory. The simulation 
of the time evolution of the chain that is large enough 
to resemble the properties in the thermodynamic limit 
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FIG. 1: The transitioning of the magnetization profile from 
the ground to the first excited state using a Not gate. Simu- 
lation obtained for a chain of 50 sites using DMRG. 



is carried out using the Density Matrix Renormalization 
Group (DMRG) algorithm. Figure [1] shows the transition 
of the magnetic profiles in the z-direction from the ground 
to the first excited state using the Not gate constructed 
from a spin-| XXZ spin chain of length 50 sites. We 
also demonstrate the construction of Pi-8, Hadamard, 
and Phase gates that form a set of universal single qubit 
gates. 

In order to have a viable quantum computing scheme 
one needs to implement at least one 2-qubit gate. Here we 
have implemented the C-Not gate which, in combination 
with the 1-qubit gates, is known to be universal Q. 

Our scheme capitalizes on the kink nature of the exci- 
tations of the XXZ Hamiltonian, which are rather sharply 
localized. We imagine a setup with two parallel chains 
with the location of the kink lined up in their ground 
states. The subspace for the 2-qubit state space is then 
I?i ® I>2 , where Vi represents the space of isolated eigen- 
values of the first chain and 2?2 for the second chain. A 
set of three controls localized near the kinks is used to 
generate the single qubit gates acting on Pi and X'2 and 
a C-Not gate on Vi (g) I?2- This scheme produces a uni- 
versal set of gates necessary for two-qubit computation. 
It is clear how to generalize this scheme to implement n- 
qubit computation. Since a universal set of single qubit 
gates and nearest neighbor C-Not gates are universal for 
n-qubit computation, this can be achieved by using n 
parallel chains and controls that are localized and act on 
neighboring chains only. 

In the next section we describe the model and review 
some of the past results. Then, in section IIIH the op- 
timal control problem to construct the quantum gates 
is described. Section |TV] is devoted to the DMRG al- 
gorithm and the specific adaptations to the XXZ spin 
chain. Finally, in section |V] we present our results based 
on numerical simulations of the XXZ Hamiltonian using 
the DMRG algorithm. 



II. THE MODEL 

In this section we describe in detail the spin-J ferro- 
magnetic XXZ model with kink boundary conditions on 
the one-dimensional lattice Z. The local Hilbert space 
for a single site a is Ha = C^-'+^ with J G = 
{0, ^, 1, |, 2, . . . }. We consider the Hilbert space for a 
finite chain on the sites [—L + 1,L] = {—L + 1, — L -|- 

2,...,+L}. This is H[_L + i,L] = The 

Hamiltonian of the spin-J XXZ model is given by equa- 
tion ([T|). Note that, by a telescoping sum, we can absorb 
the boundary fields into the local interactions: 

a=-L+l 

+ jVi - A-2 (Si - Sl+,) 

The main parameter of the model is the anisotropy A > 1 
and we get the Ising limit as A ^ 00. It is mathe- 
matically more convenient to work with the parameter 
A~^, which we then assume is in the interval [0,1]. As 
we said, A^^ — is the Ising limit, and A^^ — 1 is 
the isotropic XXX Heisenberg model. The Hamiltonian 
commutes with the total magnetization 

L 

^tot — ■ 

a=-L 

As indicated in the introduction, for each M G 
{—2JL, —2JL+ 1, . . . , 2JL}, the corresponding sector is 
defined to be the eigenspace of Sf^^^ with eigenvalue M; 
clearly, these are invariant subspaces for all the Hamil- 
tonians introduced above. These subspaces are called 
"sectors" . 

It was shown [l^, [HI, Ell that the kink boundary 
conditions lead to a family of ground states. It was also 
shown in [l^, [ll|, [l^, [l^ that for each sector there is 
a unique ground state of H^{A~^) with eigenvalue 0. 
Moreover, this ground state, ipM, is given by the following 
expression: 

«e[-L+i,L] ^ 

where the sum is over all configurations for which 
X^ct ~ ^^'^ the relationship between A > 1 and 
q G (0, 1) is given by A = {q + q^^)/2. A straightforward 
calculation shows a sharp transition in the magnetization 
from fully polarized down at the left to fully polarized up 
at the right. For this reason they are called kink ground 
states. In Q, Koma, Nachtergaele, and Starr showed 
that there is a spectral gap above each of the ground 
states in this model for all values of J. Recently Q we 
were able to prove the following theorem. 
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Theorem 1. For spin values J > 3/2, there exists a fi- 
nite Aq so that for all A > Aq, the first few excitations of 
H^{A~^) when restricted to any sector of magnetization, 
are isolated eigenvalues that persist in the thermodynamic 
limit. 

In this paper it was also proved that in eertain values 
of spin and sector, for example ^ = f and M = both 
the ground and excited states are non-degenerate (simple 
eigenvalues). This is the qubit space we work with and 
our quantum gates will be unitaries on this space. 

III. QUANTUM GATES USING QUANTUM 
CONTROL 

The problem of constructing quantum gates can be 
formulated as a problem in quantum control theory . 
The goal is to steer the system using a small number 
of control parameters such that the unitary operator de- 
scribing the quantum dynamics after a finite time T, has 
maximal overlap with a desired target unitary (the gate) . 
^From a control perspective these problems reduce to 
control of bilinear systems evolving on finite dimensional 
Lie groups. This is an optimal control problem on a two- 
level system which has been studied widely with exact 
results known in some cases. For example, time optimal 
implementation of single and two qubit quantum gates 
was studied when the Lie algebra g of su{2) (sm(4)) 
can be decomposed as a Cartan pair g = k © p with k is 
the Lie subalgebra generated by a the drift Hamiltonian 
and p is the Lie sub algebra generated by the control 
Hamiltonian's. Finding the time optimal trajectories is 
reduced to finding geodesies on the coset space G/K [G 
and K being the Lie Groups corresponding to g and k). 
The problem of driving the evolution operator while mini- 
mizing an energy- type quadratic cost was studied in p7| . 
In this case the optimal solutions can be expressed as 
Elliptic functions. The time optimal problem of popula- 
tion transfer problem of a two-level quantum system and 
bounded controls was studied in [18| and again explicit 
expressions for the optimal trajectories. In this paper we 
follow a numerical gradient based approach to optimal 
control [li,[2^. 

A. Single qubit gates 

We consider the problem of time evolution of the one- 
dimensional XXZ chain under external controls. The 
equation of motion for the unitary evolution of the 
XXZ chain isolated from the environment is given by 
Schrodingcr's equation 

U{t) = -i(Hl{A-^)+v{t)H''''''^U, f/(0) = I (2) 

In control terminology i?^(A^^) is the free or drift 
Hamiltonian and H'^^^ is the control Hamiltonian cor- 
responding to the control field v{t). We require that V is 
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FIG. 2: ErrProb, defined in Q, as a function of A ^, calcu- 
lated for the ground and first excited state i/)o and %p\, with 
^ext ^ slsl. This quantity provides a measure of the escape 
rate out of the qubit subspace. 

an invariant subspace of H'^^^ , so that the time evolution 
of the system [5] given by the unitary U{t) starting from 
an initial state in V will be constrained to V at all future 
times. The induced evolution on T) at any specified final 
time T will be the quantum gate on the qubit space T) 
and is given by the 2x2 matrix 

{U^.z)^r■^{■^^\U{TMJ) *-0,l (3) 

The control Hamiltonian we choose is the two site oper- 
ator H°^'^ = SqSi- In practice for S^Sf there is a very 
small error probability for states to move out of 1) and the 
matrix Uxxz is not exactly unitary. The matrix elements 
{■ipo\H°^^\i'k) and {iJji\H'^^^\4'k) k ^ 0,1 are proportional 
to the transition probabilities to move from states 
and ipi to other eigenstates of H^{A^^). We calculate 
the error probability to move out of the subspace V by 
the following estimates of these matrix elements 

ErrProb = II H-V^f - |(V^o|i?'="V.)P " |(V^i|i?'=''V.)l' 

(4) 

for i = 0,1. Figure [2] shows that the probabilities of 
transitioning out of the subspace V are extremely small 
for A-i < 0.3. 



B. Implementing two-qubit gates 

The idea for implementing two-qubit gates is to use 
two copies of the XXZ chain. The Hilbert space for two- 
qubit quantum computation is Pi (g) I?2 = C*, where Vi 
and P2 are the subspaces spanned by the ground state 
and first excited state of the first chain and second chain 
respectively. The Hamiltonian of an uncoupled two chain 
system is given by 

^k(^-l)(l,2) ^ ^k(^-l)(2) 

Here the notation H^{A~^Y^'> is to be interpreted as 
g^(A~^)®( I® ■ ■ ■ ® 1 ) and ff^(A-i)(2) is to be inter- 

chainl chain2 
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FIG. 3: Configuration of two XXZ chains showing the local- 
ized controls required to implement the C-Not gate. 



preted as (I ( 



The two-qubit space 



chainl chain2 

is spanned by the four vectors ipmn ■= ipm O "0" for 
771, 71 = 0, 1 which are eigenvectors of the above Hamilto- 
nian. If we consider the control system 



(5) 



U = -^[Hl{A~T^'^+v,{t){SlSf)(''> 



with U(0) = I, then by selectively turning on Vi(t) and 
V2{t) for certain time periods, the above system is equiv- 
alent to the control system ^ on chains 1 and 2 respec- 
tively during those time intervals. This can be used to 
generate single qubit gates on T>i and 1^2 ■ Moreover by 
simultaneously using vi{t) and V2{t) the local gates i.e. 
gates of the kind Xi Y2 can be generated on Vi ® I?2- 
To implement a two-qubit quantum computing scheme 
we need to also implement perfectly entangling gates i.e. 
a gate that can take a product state to a maximally en- 
tangled state. It is known that single qubit gates and 
any perfectly entangling gate are universal for two-qubit 



quantum computing 21[. Clearly such a gate cannot be 
implemented by the control scheme ([5]) alone. In this 
paper we choose to implement the C-Not gate, which 
is an example of a perfectly entangling gate. For this 
purpose we make use of an additional control namely 

We demonstrate the C-Not gate to high precision by 
using following control system 



{7 = -i(i/^(A-i)(i'2 

+V2{t){SiSf)^'^ 



(6) 



with C/(0) = I by selectively turning on and off some or 
all of the control fields vi{t), V2{t) and W3(t) for specified 



time periods. Figure |3| shows a diagrammatic represen- 
tation of the two-qubit scheme. The C-Not gate is then 
given by the 4x4 matrix with elements 



(C - Not"""), 



C. Optimal control 

We first solve the control problems ^ and ^ for the 
projected system on T> for the single chain and T>i (g) V2 
for two chain system. 



U{t) 



H + Y,Vkit)Bk)u, C/(0) 



(7) 



For the projected system on V the H and Bk's are given 
by the 2x2 matrices 



(8) 



whereas the the projected system on Vi (g) T>2 the control 
problem involves 4x4 matrices 

H.mn:rs - (tA™™ | i/^ ( ) (^^'^ |7/>,,} (9) 
{B2)mn,rs = ('/'mn I (5'o .S'l ) | T^rs ) 

where m,n^r^s = 0,1. The overlap between a desired 
unitary gate Uf and the solution of © at time T, U{T), 
is measured as the difference in the norm square ||?7/ — 
t/(r) IP, and the norm is defined in terms of the standard 
inner product :~ Tr{V'^W). The norm can be 

written as 

\\Uf - U{T)r = WUfW - 2Re{Uj\U{T)) + \\U{T)f 

and hence minimizing this norm is equivalent to maxi- 
mizing 

<P:^Re{Uf\U{T))^Tr{u}u{T)) (10) 
We define the gate fidelity as 

\Tr{U}U{T))\ 



Gate 



Tr(I) 



(11) 



To select the optimal control fields Vi(t) we use the 
numerical gradient ascent approach described in many 
books on control theory. This approach was applied to 
the quantum setting in po| . We start with the necessary 
conditions for optimality called the Pontryagin maximum 
principle which is a generalization of the Euler-Lagrange 
equations from calculus of variations. In the problems 
with costs of type (jlOp and no a priori bound on controls, 
Pontryagin's maximum principle takes the following form 
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Theorem 2. (Pontryagin maximum principle fSui. \2Bl) 
If Vi{t) 's are optimal controls of the system ^ and U(t) 
the corresponding trajectory solution, then there exists a 
nonzero operator valued Lagrange multiplier X which is 
the solution of the adjoint equations 



X{t) = -iHlt)\{t) 

a$(r) 



with terminal condition 



X{T) 



dU{T) 



and a scalar valued Hamiltonian function 
h{U{t),v,{t)) Re Tr{-iX' {t)H{t)U{t)) such that, for 
every r G (0, T] we have 



dh{U) 

dvi 



Im Tr{X {t)B,U(t)) = 



(12) 



The algorithm to find the optimal controls is as follows 

1. A suitable gate time T is chosen and discretized 
in TV equal steps of duration At = jj-. The initial 

control vf'^tk) for all the discretized time intervals 
is based on a guess or at random. 

2. For these piecewise constant controls, from C/(0) = 
I and A(r) = —[//, compute the forward and back- 
ward propagation respectively as follows 

C/M(tfc) = (13) 
AW(tfc) = ^^W(ife)FM(tfc+i)...FW(tjv)A(r)(14) 

for alHi, . . . , ijv and where r is an iteration number 
of the algorithm initially set to and 

F<^^\tk) = exp[- zAt(^H + Y,v'i"\tk)B,)} 



3. Substitute the equations and into equa- 
tion to evaluate the gradient, and then update 
the controls as 



,(^+i)(i^)=,M(,^)+,iMC/(i/c),«.(t.)) 



where t is a small step size. 



4. if ^Gate < 7 (7 being the level of accuracy) then 
done, otherwise goto step (2) for the next iteration 
with the updated controls. 

Having solved the control problem on the projected sys- 
tems to get the optimal controls vi(t), V2(t) and v^lt) we 
would like to apply them to a large system and see their 
effects on the projected system. However simulating even 
a moderately sized spin chain is hard because of the ex- 
ponentially growing dimension of the Hilbcrt space. In 
the next section we describe an algorithm by which we 
are able to simulate the XXZ chain of 50 sites. 



IV. DMRG SIMULATIONS FOR QUANTUM 
GATES 



To see the effect of the evolution of the XXZ chain 
with external magnetic controls we numerically simulate 
the XXZ chain using the DMRG algorithm. The dynam- 
ics of the interfaces of the XXZ chain using DMRG was 
studied recently in [2^ . The standard DMRG algorithm 
is a numerical algorithm originally developed by Steven 
White [13 that has worked successfully in providing very 
accurate results for ground state energies and correla- 
tion functions in strongl y c orrelated systems. Modifica- 
tions to this method [2^, [20] allow to address the physics 
of time-dependent and out of equilibrium systems. The 
crux of the DMRG algorithm is a decimation procedure 
that chooses the physically most relevant states to de- 
scribe the target states. It is now known that DMRG 
works well because the ground states of non-critical quan- 
tum chains like the XXZ chain arc only slightly entan- 
gled, i.e. they obey an area law of entanglement that says 
that the entanglement between a distinguished block of 
the chain and the rest of the chain is bounded by the 
boundary area of the block. In fact the DMRG proce- 
dure is a variational ansatz over states known as Matrix 
product states (MPS) [13 . The standard DMRG proce- 
dure and its connection with MPS and entanglement is 
described in detail in [1^ . For a single XXZ chain our tar- 
get states are the ground state V'o and first excited state 
-01 restricted to a sector of magnetization. We use the 
standard DMRG procedure with the adaptation that we 
grow the chain while restricting the blocks to the sector 
of zero magnetization using the symmetry of the Hamil- 
tonian (see [2^1). 

For the two-qubit gates we convert the two chain sys- 
tem to a one dimensional spin chain by a spin ladder 
construction. 



n 



(1) 

~L+1 



n 



(1) 

-L+2 



H 



(1) 



H 



(1) 

[-L+1,L] 



0,(2) _ ^(2) 

'T-L — 'T-f-L+l.Ll 



The single site Hilbert space for the DMRG is the rung 

composed of Ha'' ®'Ha'' for a G [— L + 1..L]. On this site 
we define the local operators 

ST^=Sl®la, ST^ = la®Sl for i= 1,2,3 

We can then write the Hamiltonian of this single chain 
using the above construction 

Hl{A-^)^''^-^ := X: <Ui(A-^) + /^i^,Ui(A-i) 

a=-L+l 

for k = 1,2. We carry out the DMRG procedure 
as described in the algorithm with the Hamiltonian 
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H^{A~^Y^''^^ but we ensure that we keep both the chains 
in the magnetization sector by simultaneously diagonal- 
izing H^{A~^)'--^'^'> with the total magnetization opera- 
tors 

E Sl^'^^ fork=l,2 

a=-L+l 

The target states ipo^ipo, ipo'^ipi, V'l ® V'o , "01® ipi are 
the simultaneous eigenvectors of the these operators and 
form the computational basis |00), |01),|10) and |11) for 
two-qubit quantum computation. 

To compute the time evolution of the chain under the 
controlled evolution by the external fields we use the 
time-dependent DMRG procedure. The idea is that a 
two site operator can be applied to a DMRG state most 
effectively by expressing the state in the basis where the 
left block has length a; — 1 so the two middle sites that 
are untruncated are the the sites where the operator is 
acting. We can write the time evolution in the Trotter 
decomposition 

. . . g-^/ii,-2,t-ig-ifti.-i.i, 0(^3^ 

To apply e^'^* to the ground and excited states in the 
basis with the center sites all the way to the left we ap- 
ply e~2''-i'+i.-J^+2. After shifting one site to the right we 
apply e~5''-i'+2,-i.+3 gtc. Since all our controls are two 
site controls at the center, only the interaction /io,i is 
time-dependent. In the adaptive time-dependent meth- 
ods the Hilbert space is continuously modified as time 
progresses by carrying out reduced basis transformations 
on the evolved state. In our case since the gates are ob- 
tained in a relatively short period of time our Hilbert 
space remains unchanged resembling the static DMRG 
methods. 



V. RESULTS 

In this section we present numerical results of the con- 
struction of quantum gates using the spin-3/2 XXZ spin 
chain. Our results are for the universal set of single qubit 
gates consisting of the Not (X), Hadamard (H), Pi-8 (T) 
and Phase (S) gates and the two-qubit C-Not gate. All 
results are obtained using the DMRG algorithm and the 
optimal control methods described in the previous sec- 
tions. 

The steps carried out to obtain the single qubit gates 
are as follows 

1. We use ground state DMRG of the XXZ chain 
to obtain the lowest eigenvectors ■00 and ipi of 
H^{A^^) in the sector corresponding to M = 0. 

2. We obtain the projected 2x2 control system of 
equation ^ with matrices H and Bi with matrix 



IN Ot yJ\. J 
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U.4o00 


-u.Doyy 


A Q /I /I 1 

-U.o44i 


A OAdd 

-U.z4D0 


u.zooy 


-U.Dz4i 


A QOOO 


A 01 AQ 

-U.ziUo 


n AO A fi 

-U.Uz4d 


-u.ouyy 


A oc:nA 

-U.zoyu 


A 1 CQQ 

-U.iooo 


-U.zDSl 


-U.o ( zo 


A 1 QOQ 

-U.lozo 


A 071 K 

-U.U / lo 


n A QOO 

-U.4oyy 


-U.z4U4 


u.u ( uy 


U.UooD 


-0.5065 


-0.1424 


0.3368 


0.1704 


-0.4553 


-0.1001 


0.5877 


0.3053 


-0.2959 


-0.1225 


0.7029 


0.4128 


-0.0605 


-0.2033 


0.6294 


0.4642 


0.1909 


-0.3236 


0.4374 


0.4516 



TABLE I: Numerical simulation of the construction of the 
Not, Pi-8, Hadamard and Phase gates. Results are obtained 
using DMRG and time-dependent DMRG for a spin-| chain 
of L = 50 sites and at = 0.3 in the sector corresponding 
to M = 0. The table shows the values of the control field 
vi{t) with gate time T = 10 discretized with At = 0.5. 




FIG. 4: (Clockwise) The controls for Not, Hadamard, Pi-8 
and Phase gates controls plotted versus time discretized for 
time steps of At = 0.5. 
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V2(t) 




-0.4040 


0.3953 


0.0540 


2.4494 


0.2588 


0.0501 


3.7163 


0.1886 


-0.1314 


3.0455 


0.1766 


-0.2677 


1.6565 


0.2185 


-0.0872 


0.5583 


0.3085 


0.1455 


-0.1036 


0.4206 


0.2346 


-0.5117 


0.5196 


0.1648 


-0.8215 


0.5716 


-0.0282 


-1.0083 


0.5578 


-0.3599 


-0.8630 


0.4926 


-0.9102 


-0.3774 


0.3612 


-1.4955 


-0.0810 


0.0978 


-1.4905 


-0.0940 


-0.1214 


-0.9196 


-0.0881 


-0.1138 


-0.3217 


-0.0408 


0.0794 


0.0922 


-0.1805 


0.3399 


0.3374 


-0.8094 


0.5868 


0.4146 


-2.2183 


0.8522 


0.2021 


-4.2425 


1.2694 


-0.4760 



TABLE II: The C-Not gate controls using two spin-| XXZ- 
chains of length L = 50 at = 0.25 and M = for both 
the chains. The gate time T = 3.5 is discretized into A*' = 20 
time steps. The table shows the values for the three control 
fields vi , V2 and vs that are constant during any one of the 
time intervals. 



The gates obtained using the XXZ chain and their fideh- 
tics are as follows. 



0.0016 -O.OOlli 0.0033 + 0.9997* 
^-0.0017 + 0.9997i 0.0017 + O.OOlli^ 
J'x = 0.9997, 



-0.0027 + 0.7081i 0.0011 + 0.7053i 
^-0.0016 + 0.7052i -0.0022 - 0.7085i^ 
Th = 0.9995, 



0.9221 - 0.3859i -0.0037 -f 0.0038i 
, 0.0037 -F0.0038i 0.9216 0.3871i , 



= 0.9995, 



Sxxz ~ 

Ts = 0.9997 



0.7043 - 0.7095i -0.0046 + 0.0015i 
0.0045 + 0.0016i 0.7017 + 0.712H , 



elements given by ([5]). For a target gate Uf and 
a suitable final time T we find the optimal control 
Vi{t) on this 2x2 system using the technique de- 
scribed in section UlI CI 

3. Finally we apply the time-dependent DMRG pro- 
cedure of section HVl to the chain of ^ for a speci- 
fied time T starting from iJjq and ipi and using the 
vi {t) found in step 2 to get the time evolved states 
^o(T) = f/(T)^/>o and = U{T)i'i. We com- 

pute the induced evolution on the subspacc V to 
obtain the gate Uxxz given by the matrix elements 
{ipi\tljj{T)) for i,j = 0,1 and compare the overlap 
with Uf using equation (jlip . 

Our desired single qubit target gates are given by the 
unitarics 




The optimal controls vi {t) used to get the gate results are 
shown in Table [J and Figure HI For the C-Not gate the 
procedure described earlier is only slightly modified. We 
do the ground state DMRG of a one dimensional chain 
built from the spin ladder described in section IIVI to get 
four eigenvectors ipmn for m, n = 0, 1. The optimal con- 
trol procedure is applied to the 4x4 control system ([7]) 
with H, Bi, i?2, B3 given by equations Q to find the con- 
trols vi{t), V2{t) and V3{t). The time-dependent DMRG 
procedure is applied to the chain of equation ([6]) for time 
T with the controls vi{t), V2(t) and v^lt) to get the time 
evolved states ipmniT) ~ U{T)ijjmn- The induced evolu- 
tion on the subspace Vi (g) V2 gives the C — Notj;^;^ gate 
with matrix elements {'4'mn\'4^rs{T)). Table Ull shows the 
optimal controls vi{t), V2{t) and V3{t) used to obtain the 
C-not gate. The gate obtained using the XXZ chain and 
gate fidelity is as follows 



8 



C-Not 



C - Not^ 



/l o\ 

10 

1 
\0 1 0/ 

/ 0.9959 + 0.0001i 



-0.0015 + 0.0006i 0.0003 - O.OOOSi -0.0010 - O.OOOlA 



-Not 



-0.0014 -O.OOlOi 0.9939 - O.OOOSi 0.0015 + O.OOOOi -0.0005 + 0.0005i 

0.0013 -O.OOOH 0.0004 + 0.0003i 0.0004 - 0.0003i 0.9945 - O.OOOSi 

y-0.0003 - 0.0003i -0.0013 + 0.0002i 0.9954 - 0.0004i 0.0003 + 0.0003i / 
0.9949 
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